APPROXIMATIONS  TO  THE  DIRECTIVITY 
INDEX 


Robert  J.  Bobber 

Naval  Research  Laboratory 
Washington,  D.  C. 

1  May  1974 


DiSTftIBUTED  BY: 


mntiii  iiciwcii  miiMiiii  Mificc 

U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road.  Spiingfieid  Va.  22151 


i 

t: 

i 

r 


fe 

t 


j 

i- 


fc." 


St 


- WCLftSSiriED -  AK'7'7(////>> 

SECORITY  classification  of  this  face  nwi»n  OMa  EfitnaiO  H"  D  /  f  fj  L5»^C/ 

REPPUT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  NEPONT  number  2.  govt  accession  NO. 

NRL  Report  7'’50 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TlTLt  (md  SubttiU) 

APPROXIMATIONS  TO  THE  DIRECTIVITY  INDEX 

S.  type  of  report  a  PERIOO  COVERED 

Final  report  on  one  peu±  of 
the  problem 

S.  performing  ORG.  report  NUMBER 

7.  authori*; 

Robert  J.  Bobber 

t.  CONTHACT  ON  GRANT  NUNBENf.J 

3.  performing  organization  name  ano  aooress 

Naval  Research  Laboratory 

Underwater  Sound  Reference  Division 

J.O.  Box  8337,  Orlando-  FL  32806 

10.  PROGNAN  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

NRL  Problem  K03-30.401 

It.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Depeurtnent  of  the  Navy 

Office  of  Naval  Research 

Arlinoton.  VA  22217 

13.  REPORT  DATE 

1  May  1974 

IS.  number  OF  PACES 

iii  +  17 

t4.  monitoring  agency  Ni^ME  A  AOORESVIl  dtHdtmt  tnm  Contnlthig  Otitem) 

IS.  SECURITY  CLASS,  (ot  Ala  rapertj 

UNCIASSIFIED 

1S.L  OECLASSiFiCATtON/DONNCRAOINC 
SCHEDULE 

IS.  CISTRIBUTION  statement  (oI  Ala  RaRoMJ 

Approved  for  public  release;  distribution  unlimited 

17.  OlsTRIBUTION  STATEMENT  (ot  Aa  aSatracI  aniarad  m  Slack  39,  1/  dllfanni  turn  SaFactJ 

IS.  supplementary  notes 

fi 

D  D  C 

nPjn^ 

If.  KEY  WORDS  fCoiitMiM  mn  mfdd  it  ndcmmddfy  mtd  Idmtfty  kr  btmck  mm^m)  j 

_ .  ^  ,  Rcprosluccd  by 

Directxvity  index  NATIONAL  TECHNICAL 

U.iderwater  sound  transducers  INFORMATION  SERVICE 

0  S  Department  of  Commerve 

SprtnsfienJ  VA  2215* 

:[<  MAY  15  1.974 

llEIEtSrO  LSlUJ 

D 

20.  ABSTRACT  {Cdtitimf  mn  frdt—  mtdd  it  **4  ithntiff  ky  mmtbde) 

Directivity  indexes  of  three  different  underwater  sound  transducers  have 
been  measured  or  computed  by  five  different  methods;  (1)  theoretical  calcula¬ 
tions  using  piston  area  or  line  length,  (2)  beam-width  measurement  and  theo¬ 
retical  calculation,  (3)  pattern  raessurement  and  graphical  integration,  (4) 
direct  digital  directivity  ii.dex  measuring  3y.stera,  and  (5)  diffuse-.sound  method 
The  resiilts  indicate  that  the  directivity  index  of  emy  ordinary  transducer  can 
be  obtained  fr«n  calculations  based  on  known  configurations  and  dimensions  or 

DO 


f  SUM 
I  JAN  73 


:m 


COITION  OF  I  NOV  Cl  IS  OBSOLCTC 
S/N  0103>01€>eS0t  I 


UlfCLASSIFIEO 


i 


SCCuntTr  classification  of  TKIS  FAOC  (NNm  Oala  CMMO 


_ UNCLASSIFIED _ 

.uuumTV  CUASSIrlCATIOM  or  this  PACEf»/i«ii  Df  gwtfO 


20.  Abstract  (continued) 

.  beam-width  measurements  with  a  degree  of  reliability  and  accuracy  that  is  no 
worse  than  any  measurement  technique  and  that,  in  roost  cases,  elaborate 
roeasurements  for  determining  directivity  index  aire  not  justified. 


_ UNCLASSIFIED _ 

security  CLAMiriCATtON  OF  TKIS  RACErVMn  DM*  JblMO 


Contents 


Introduction  .  i 

Comparison  of  Five  Methods  .  .  i 

Energy  in  Minor  Lobes  .  4 

Circular  Pistons  .  5 

Rectangular  and  Other  Pistons  .  7 

Radiating  Areas  .  8 

Beam  Widths  .  9 

Rear  Lobes  .  9 

Radiation  Resistance  .  lo 

Line  Sources  .  ii 

Shaded  Transducers  .  13 

Arrays  with  Mutual  Coupling  .  15 

Conclusion . 16 

Acknowledgement . 16 

References . 16 

Figures 

1.  Array  configurations  for  (a)  one  quadrant  of  typs  F27 
transducer,  (b)  one  quadrant  of  type  F33  transducer,  ^md 

(c)  type  F37  line  transducer .  2 

2.  The  directivity  index  of  (a)  type  F27  transducer,  (b) 
type  F33  transducer,  and  (c)  type  F37  transducer,  as 

determined  by  methods  1,  2,  3,  4,  and  5 .  4 

3.  The  directivity  index  of  a  circulsu:  piston  of  diameter  d 

without  baffle,  with  a  rigid  plane  baffle,  in  the  end  of 
an  infinitely  long  rigid  pipe,  and  as  approximated  by  the 
expression  10  log  (4irA/x2) .  6 

4.  The  directivity  index  of  a  rectangular  piston  in  an 
infinite  plane  rigid  baffle  as  a  fraction  of  length  L 

or  width  W .  7 

5.  The  directivity  index  of  a  line  source  according  to 
Stenzel's  exact  expression  (Eq.  (8)),  and  three  degrees 

of  approximation  by  Eq.  (10) . 12 

6.  Major  lobes  of  patterns  for  a  uniform  line,  a  linearly 

tapered  line,  and  10-elcment  binoaially  shaded  line  array  ....  15 

Tables 

1 .  Percentage  of  acoustic  energy  in  each  lobe  of  an  idead 
piston  with  a  diameter  of  d/X  %ravelengths ,  and  an  ideal 

line  with  a  length  of  L/X  wavelengths .  5 

2.  Correction  (In  decibels)  to  calculated  directivity  index 

as  a  functira  of  rear  lobe  height . 10 


iii 


APPROXIMATIONS  TO  THE  DIRECTIVITY  INDEX 


Introduction 

The  directivity  index  has  long  been  an  ii!5>ortant  parameter  in  eval¬ 
uating  some  types  of  electroacoustic  tramsducers — particulaurly  in  sonar 
applications.  It  is  a  difficult  parameter  to  analyze  and  measure  because 
a  three-dimens:Lonal  integration  of  the  radiation  pattern  is  involved. 

For  this  reason,  approximations  have  been  widely  used  in  both  the  theory 
and  measurement.  In  spite  of  the  fact  that  Stenzel  [1] ,  Holloy  [2] ,  and 
others  worked  out  the  fundamental  theory  many  years  ago,  there  hcis  been 
little  quemtitative  analysis  of  the  theoretical,  eind  measurement  apprcxi- 
mations.  At  the  same  time,  trainsducer  designers  have  either  gone  to 
elaborate  and  costly  techniques  to  measure  the  directivity  index  [3-8] , 
or  have  used  *ry  single  ccxnputat.'onal  aids  such  as  special  slide  rules 
that  are  based  on  idealized  models  [9] .  The  purpose  of  this  study  was 
to  cooq>are  the  results  of  these  two  extremes  in  methodology,  to  quamtify 
some  of  the  limits  of  approximations,  to  identify  the  most  feasible 
method  of  detenuining  the  directivity  index,  and  to  report  some  results 
of  two  little  known  or  used  measurement  methods.  An  tinexpected  side 
light  has  been  the  identification  of  am  error  in  Stenzel 's  original 
analysis  amd  scxne  errors  in  translating  Stenzel ‘s  work  into  English 
[10,11]  amd  in  preparing  a  second  edition  of  his  book  [12]. 

The  principal  conclusion  from  this  study  was  that,  in  most  cases, 
elaborate  measiirements  for  determining  the  directivity  index  are  not 
justified  by  the  accuracy  or  precision  of  the  results. 


Cm^arison  of  Five  Methods 

D.'jrectivity  factors  of  three  different  underwater  sound  transducers 
have  been  measured  or  coc^ted  by  five  different  methods. 

The  three  transducers  are  NRL-OSRD  types  F27,  F33,  and  F37.  The 
F27  apprordmates  a  uniform  circular  piston  in  a  rigid  baffle,  as  shown 
in  Fig.  la.  It  is  cocaprised  of  an  eunray  of  55  lead  metaniobate  disks. 
The  type  F33  approximates  a  nonuniform  circular  piston  in  a  rigid  baffle, 
as  shown  in  Fig.  lb.  It  is  coBq>rised  of  two  arrays.  The  outer  array 
contains  64  barium  titanate  rectangiilar  plates;  the  inner  2unray,  12  lead 


The  five  xoethods  are: 

1.  Theoretical  Calculations  Using  Pist.m  Area  or  Line  Length.  This 

iDethod  requires  only  theoretical  calculatir-^  based  on  well-known  approxi- 
oations.  The  directivity  factor  of  a  piston  in  an  infinite  rigid  baffle 
is  given  approximately  by  the  expression  where  A  is  the  piston 

area  and  1  is  the  wavelength.  Foj  a  line  or  thin  cylinder,  the  corre¬ 
sponding  expression  is  2L/X,  where  L  is  the  length.  These  are  single 
expressions,  but  uncertainties  usu^d.ly  arise  in  ascertaining  A  and  L 
because,  in  practice,  both  pistons  and  lines  eure  really  eurrays  of  elements. 
The  spaces  between  the  elements  usually  are  included  in  A  and  L,  but  the 
effective  edges  or  ends  of  the  arrajs  are  sore  indefinite,  and  generally 
introduce  an  uncertainty  of  about  5%  for  L  and  10%  for  A. 

2.  Beam-Vidtk  Measurement  ar.d  Theoretical  Calculation.  The  beam  width 
of  the  r^iiation  pattern  was  measured.  The  transducer  size  <tnd  shape  were 
inferred  from  these  seetsureoents ,  and  the  directivity  factor  then  was  ceQ.- 
culated  in  a  manner  3imil^ur  to  the  first  method.  In  both  calculation 
methods,  it  is  eissumed  that  the  differences  between  the  theoretical  and 
rezd.  minor  lobe  structures  in  the  pattern  are  negligible. 


^Wjl!l|lPW''^i(iW*i' ,  iwif  »<“!•,'  ">1'  ‘  '«|»W  Wi(|()Mtffi.|ii')qip;i)||i{)||)||i|M  •f.-T  ii|,t*rt*W||WWV''‘' 


3.  Pattezm  Measurement  and  Graphical  Integration.  A  number  of  two- 
dimensional  patterns  were  plotted.  Then  graphical  integrations  were  carried 
out  to  ascertain  aK>roxiinately  the  three-dimensioncil  pattern.  This  is  a 
stamdcird  metliod,  but  very  time  consuming  xinless  the  pattern  has  circulcir 
symmetry  about  at  least  one  axis,  or  some  con^terized  technique  [4,5]  is 
used  for  the  integration. 

4.  Direct  Digital  Directivity  Index  Measicring  System.  A  new  digital 

»  system  was  used  [6,7] .  A  seven-element  semicirculcir  hydrophone  array  was 

swept  through  a  spherical  surface  around  the  tramsducer,  and  25;'  sarnies 
of  the  raidiated  sound  pressure  level  were  obtained  in  a  few  minutes.  The 
252  values  were  processed  by  digital  consular  methods  to  obtain  the  direc¬ 
tivity  factor  in  a  short  time.  Wo  patterns,  per  se ,  are  required  for  this 
method. 

5.  Diffuse-Sound  Method.  The  identity  between  the  directivity  factor 
amd  the  ratio  of  the  free-field  to  the  difiuse-field  receiving  sensitivity 
of  a  tramsducer  [13]  was  used.  Diffuse  fields  or  reverberent  chambers 
have  been  little  used  in  underwater  acoustics  because  of  the  long  wave¬ 
lengths  and  difficulties  in  obtaining  laurge  iii^>edances  mismatches.  The 
diffuse  field  sensitivities  tised  in  this  experiioent  were  obtained  by  B.  G. 
Watters  in  the  reverberation  tank  at  the  Bolt,  Beranek,  and  He%iman  coii^>amy 
i-  Catabridge,  Haiss.  [14].  This  tank  has  the  dimensions  9.75  x  7.01  x  4.27 
me't-'rs  and  reverberation  times  an  long  as  5  seconds.  The  only  other  re¬ 
po;,  ted  use  of  this  method  is  from  Reznikov  amd  Snytko  [15]  who  used  both 
spatial  and  tenoral  averaging  in  a  smatLl  water-filled  vessel.  Their  work 
is  difficult  to  assess  or  \ise  because  of  an  incoc^lete  description  of  the 
transducer,  the  use  of  unexplained  corrections,  and  results  which  show 
that  the  directivity  factor  of  a  cylindrical  transducer  is  not  proportional 
to  frequency,  as  it  should  be. 

The  directivity  index,  or  ten  times  the  logaurithm  of  the  directivity 
factor,  is  shown  in  Fig.  2  for  the  three  transducers  as  determined  by  the 
five  methods.  The  calculated  directivity  indexes  for  the  types  F27  amd 
F33  from  the  first  method  au:e  shown  am  broad  lines,  0.5  dB  wide,  because 
of  the  uncertainty  in  the  vad.ue  to  be  used  for  the  aurea  A. 

Che  a'  reement  among  all  methods  except  the  fifth  (diffuse  field)  is 
imus^uall'.  good  for  the  F33 — so  good  that  all  data  points  fall  within  the 
0.5-dB  spread  of  the  calciilated  vad.ues.  The  diffuse  field  data  are  clearly 
too  liigh.  The  discxepamcy  probably  is  due  to  the  is^erfecc  diffuseness  of 
the  field. 

The  scatter  among  the  methods  is  greater  for  the  F27,  but  the  average 
of  the  three  experimental,  methods  (3,  4,  and  5}  agrees  well  with  the  two 
cailculation  methods  (1  amd  2). 

For  the  F37,  four  of  the  five  methods  aure  in  good  agreement  above 
20  'kHz,,  but  this  time  it  is  the  digital  method  that  does  not  agree;  how¬ 
ever,  cnly  one  (25  kHz)  of  the  three  data  points  is  widely  different. 

Below  20  )cHz  the  diffuse  field  method  is  agadn  too  high. 
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10  IS  20 

Fr«<}uency  (kHz) 

(b) 

Fig.  2.  The  directivity  index  of  (a)  type  F27  transducer, 

(b)  type  F33  transducer,  and  (c)  type  F37  tramsducer,  as  de¬ 
termined  by  methods  1,  2,  3,  4,  a.ad  5. 

The  data  show  that  the  diffuse-field  method  in  the  BBN  tank  should  be 
liaj.ted  to  25  kHz  and  higher  frequencies.  Eeyond  that,  there  is  not  suf¬ 
ficient  consistency  in  the  data  to  conclude  that  one  of  the  three  experi¬ 
mental  methods  is  to  be  preferred. 

llie  one  conclusitm  that  emerges  from  the  ea^riment  is  that  the  direc¬ 
tivity  index  obtained  from  the  theoretical  or  from  the  beam  width  calcu¬ 
lations  is  as  reliable  as  any  of  the  experimental  methods,  and  it  appears 
futile  to  go  to  elaborate  measurements  for  transducers  of  conventional 
shapes . 


Energy  in  Minor  Lobes 

Table  1  shows  the  acousti*:  energy  distribution  among  the  major  cuid 
minor  lobes  in  typical  patterns  It  is  evident  from  this  table  that  minor 
lobes  contribute  very  little  to  the  directivity  factor  or  inde:s.  Neglect¬ 
ing  all  the  minor  lobes  would  introduce  an  error  of  less  than  0.8  dB  for 
a  pistem  and  less  than  0.4  dB  xor  a  line.  In  practice,  of  course,  it  is 
not  a  matter  of  entirely  neglecting  the  minor  Ic^s,  but  rather  neglecting 
the  difference  between  the  idee^-ized  and  the  real  pattern.  Clearly,  1- 
or  2-d3  variations  in  the  height  of  the  first  minor  lobe  aud  even  larger 
variations  in  the  others  are  not  going  to  make  perceptille  differences 
between  the  real  directivity  index  and  the  directivity  -ndex  based  on  an 
ideal  model  and  measurements  of  only  the  major  lobe.  This  fact,  of  course 
supports  the  thesis  that  once  the  beam  width  of  the  major  lobe  is  knewn, 
along  with  the  basic  configuration  of  the  radiator  (circle,  rectangle, 
cylinder,  etc.),  and  the  knowledge  that  the  minor  lobe  structure  is  net 
radically  abnormal,  no  further  neasiirement  is  necessary. 
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Table  1.  Percentage  of  acou<itic  energy  in  each  of  am  ideal 
piston  with  a  diaMter  of  d/A  wavelengths,  and  an  ideal  line  with  a 
length  of  L/A  wavelenv^ths. 


Minor 

Leases 


d/A 

n/A 

Major 

hotoe 

1 

3 

4 

5 

6 

7 

2 

0.85 

0.15 

3 

.84 

.09 

0.07 

5 

.84 

.08 

.03 

0.02 

0.03 

8 

.84 

.08 

.03 

.02 

.01 

0.00 

0.00 

0.02 

2 

0.95 

0.05 

3 

.93 

.05 

0.02 

5 

.92 

.05 

.02 

0.01 

0.00 

8 

.92 

.05 

.01 

O 

4 

.00 

0.00 

0.00 

0.00 

An  interesting  side  light  is  shown  in  Table  1.  The  energy  ^n  suc¬ 
cessive  side  lobes  does  not  always  diminish  steadily.  In  large  pistons 
the  last  minor  lobe  contains  more  energy  them  some  intermediate  lobes. 
The  larger  solid  angle  of  the  last  lobe  more  than  compensates  for  the 
lower  average  level. 


Circular  Pistons 

After  one  is  persuaded  that  ceilculations  based  on  approximations  are 
sufficient  in  most  cases  to  determine  the  directivity  index,  it  is  stilL 
necessary  to  define  quantiteti\'e  limits  for  these  approximations.  Of  the 
coanon  configurations,  the  circulsu:  piston  in  a  rigid  plane  baffle  is  the 
best  known.  Figure  3  shows  the  directivity  Lidex  for  the  rigorous  case 

coeputed  frea 


5^  -  10  log 


(lcd/2)2 

2JiCcd) 

kd 


where  k  *  2»/A,  d  is  the  diameter  of  the  piston,  and  denotes  the  first- 
order  Bessel  functicxi,  and  for  the  ai^roximation  based  on  the  area. 
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-  10  log  (4iiA/A^).  (2) 

It  also  shows  the  directivity  index  from  Beranek  [16]  for  the  two  most 
coBDon  departures  from  am  infinite  plane  rigid  baffle.  Many  underwater 
sotmd  transducers  approximate  the  piston  in  the  end  of  a  long  tube.  An 
tuibaffled  loudspeedcer  is  an  example  of  the  \inbaffled  piston,  which  approxi 
mates  a  dipole  with  a  directivity  index  of  4.7  dB  at  low  frequencies  and 
is  consistently  3  dB  lower  tham  the  baffled  piston  at  high  frequencies 
because  of  its  bidirectional  pattern.  At  high  frequencies,  it  is  the 
same  as  a  plane  baffled  piston  radiating  in  both  directions. 

The  directivity  index  usually  is  a  useful  patrameter  only  when  the 
major  lobe  is  sonewhat  narro:«  amd  the  index  is  of  the  order  of  10  dB  or 
more.  From  Figure  3,  Eq.  (2)  is  cleaurly  a  very  good  approximation  for 


Pipe  btffle 


Pl*ne  Uffl» 


No  baffle 


10  log  {4tA/i‘) 


Fig.  3.  The  directivity  index  of  a  circulau:  piston  of  diameter 
d  without  baiffle,  with  a  rigid  plane  baffle,  in  the  end  of  am 
infinitely  long  rigid  pipe,  and  as  ac^roximated  by  the  expression 
10  log  (4rA/l2). 
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diameters  of  several  wavelengths  with  either  type  of  baffle,  and  even  with 
no  baffle  if  the  3-dB  correction  is  subtracted.  The  lower  limit  of  the 
approximation,  for  baffled  pistons,  is  d  =  X  for  errors  less  than  0.1  dB, 
and  d/2  >  X  for  errors  less  than  1.0  dB. 


Rectangular  and  Other  Pistons 

Theoretical  values  for  the  directivity  index  of  rectangular  pistons 
in  infinite  plcme  rigid  baffles  have  been  calculated  by  both  Stenzel  [1] 
euid  Molloy  [2].  Stenzel  [17]  has  pointed  out  some  errors  in  Molloy's 
paper,  so  Stenzel 's  values  have  been  used  in  preparing  Fig.  4  whenever  the 
two  authors  disagree.  The  approximation  based  on  the  area  of  the  rectangle 


stenzel  t  Nolioy 
10  log 


Beffled  line 


Fig.  4.  The  directivity  inder  of  a  rectangular  piston  in  2m  infinite 
plane  rigid  baffle  as  a  function  of  length  L  or  width  W.  Solid  lines: 
2m  determined  frco  the  approximation  b2med  on  2u:ea  [d^  s  10  log 

(4ri;W/X2)J.  Demhed  lines:  t2d(en  from  data  by  Stenzel  or  Moiloy. 
Dzmh-dot  line:  from  Stenzel,  2ifter  applying  a  3-dB  correction  for  a 
baiffled  line  radiating  into  a  h2d.f-space. 
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is  cilso  shovm  in  Fig.  4.  From  these  data,  it  is  appeu  oit  that  the  limits 
for  the  approximation  cure  at  least  as  low  eis  for  the  c;ircular  piston.  That 
is,  where  both  L  and  W  equal  or  exceed  a  wavelength,  or  kL/X  ^  r  eind 
kW/X  ^  r,  the  approximation  error  is  negligible.  And  where  both  L  amd  W 
equad  or  exceed  a  half  wavelength,  the  error  is  less  than  1.0  dB.  It  is 
interesting  to  note  than  when  kW/2  =  1  or  W/X  =  l/ir,  the  approximation 
AnA/X^  reduces  to  4L/X,  or  the  approximation  for  a  line  source  in  a  baffle 
radiating  into  a  half-space. 

It  is  evident  from  Figs.  3  and  4  that  a  badfled  piston  of  any  shape 
intermediate  between  a  circle  and  rectangle  (ellipse,  octagon,  etc.) 
would  have  a  directivity  index  accurately  given  by  10  log  4iiA/X^,  provided 
that  its  smallest  dimension  is  one  wavelength.  It  is  also  evident  that 
for  dimensions  greater  than  a  wavelength,  the  baffle  configuration  makes 
no  difference.  Further,  the  directivity  index  of  a  piston  of  some  unusual 
shape  (a  cross,  for  exaaq>le)  also  is  given  by  10  lo^  4iiA/X^,  provided  that 
its  narro^fest  dimension  exceeds  a  wavelength.  This  result  follows  from 
the  argument  that  if  the  piston  area  is  subdivided  into  segments,  and  the 
dimensionst  of  each  individual  segment  meet  the  wavelength  criteria,  then 
the  whole  radiator  meets  the  apprrjcimation  criterion. 


Radiating  Areas 

Given  that  Eq.  (2)  is  a  valid  theoretical  approximaticax  for  rsost  pis¬ 
ton  transducers,  tliere  still  remains  the  problem  of  determining  the  piston 
area  A  in  real  transducers  as  illustrated  by  Fig.  1.  It  is  a  rule  of 
thumb  in  sonu  transducer  design  that  if  the  element  spacing  in  an  array 
does  not  exceed  0.8X,  the  array  then  functions  essentiaQ.ly  as  a  plane 
radiatoi;.  Insofeu:  as  the  directivity  pattern  is  concerned,  this  meani; 
that  the  major  1«^  of  the  array  is  the  same  as  if  the  array  were  a  uni¬ 
form  pl:me.  The  minor  lobes,  however,  are  quite  different  until  the 
spacing  becomes  less  than  0.2X  [18].  The  use  of  the  0.8X  rule  is  asply 
supported  by  the  iiiq>lication  from  Table  1  that  deviations  in  the  minor 
lobes  Ciin  be  neglected  for  purposes  of  determining  radiated  energy. 

Within  the  limit  of  the  0.8X  rule,  the  interstitial  spaces  in  an  array 
are  included  in  the  theoreticeO.  radiating  area. 

The  periphery  of  the  array  is  more  of  a  problem.  It  would  seen  logi¬ 
cal  that  heilf  an  interstitieO.  space  coepletely  surrounding  each  eles^t 
should  be  included  in  the  uea.  This  adds  a  thin  periph  ral  area  that, 
in  the  case  of  the  F33  shown  in  Fig.  lb,  is  a  uniform  thin  border  a  half- 
interstitiad-space  wide.  But  for  the  F27,  shown  in  Fig.  la,  it  is  not  so 
straightforward  because  of  the  unusual  shape  of  the  interstice.  Further, 
if  the  interstices  are  included  in  the  area,  why  not  some  of  the  concave 
comers  at  the  periphery?  Calculatitxis  of  the  effective  aureas  of  the  F27 
and  F33  were  maide  in  various  ways,  including  subjective  judgements  in 
soae  cases.  Tl»  results  showed  a  spread  of  about  10%  in  the  area,  or 
0.5  dB  in  the  directivity  index.  An  area  detexsdned  by  averaging  the 
results  of  several  techniques  is  probably  the  only  practical  method. 
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t..dm  Widths 

The  effective  area  of  a  piston  can  be  found  by  measuring  tl’^e  beam 
vjidth  of  the  major  lobe  in  the  pattern,  provided  the  area  is  or  ;pproxi- 
mates  a  circle,  squaire,  or  rectangle. 

The  diameter-to-wavelength  ratio  d/X  of  a  circular  piston  is  given  in 
terns  of  the  6-dB-down  half  beam  width  6  by 

d/X  =  0.70/sin  6.  (31 

Similarly,  the  side  of  a  square  or  rectemgle,  or  the  length  of  a  line, 
is  given  by 

L/X  =  0.60/sin  6,  (4) 

where  0  is  the  6-dB-down  half  beam  width  of  the  pattern  in  the  plame  of 
the  dimension  L. 

The  directivity  index  then  is  found  from  the  dimensions  au:d  Eq.  (2) . 

This  method  has  the  advantage  of  dealing  directly  with  the  radiated 
acoustic  energy.  The  disadvantage  is  that  most  transducers  do  not  have 
ideal  shapes,  as— for  example— the  F2?  and  F33,  and  beam  widths  in  several 
planes  must  be  averaged.  This  was  done  for  the  data  shown  for  method  2 
in  Fig.  2.  The  results;  in  Fig.  2  indicate  that  averaging  beam  widths  pro¬ 
duces  about  the  same  results  as  averaging  areas. 

Beam  widths  are  measured  at  either  the  3-,  6-,  or  10-dB-down  points. 

It  was  found  that  using  the  and  10-dB  down  values  for  0  in  Eqs.  (3) 
and  (4)  produced  the  most  consistent  results,  though  using  the  3-dB-down 
bean  width  produced  directivity  indexes  only  0.1  or  0.2  dB  different  from 
the  other  two.  It  probably  is  ,i>est  to  measure  all  three  beam  widths  and 
check  against  the  theoretical  values  that  show  ths  relative  beam  widths  as 

3  dB  down  0.73 

6  dB  down  1.00 

10  dB  down  1.23 

These  ratios  apply  to  both  circular  pistons  and  square  or  rectcuigular 
pistons  (and  lines)  in  planes  parallel  to  a  side. 

The  fact  that  the  ratios  are  the  same  for  both  illustrates  that  the 
two-dimensional  patterns  have  the  same  relative  shape  in  the  direction 
where  most  of  the  sound  energy  is  ladiated. 


Rear  Lobes 

Perhaps  the  mort  cons^.‘^^ent  difference  between  the  patterns  of  real 
imd  ideal  piston  transducers  is  in  the  existence  of  rear  lobes.  Ideally, 
there  would  be  no  rear  lobes.  In  practice,  they  appear  often  because  it 
is  so  difficult  to  obtain  a  truly  rigid  baffle  or  housing  in  underwater 
acoustics . 


If  the  reau:  or  back  plate  of  a  transducer  housing  vibrates,  a  pattern 
lobe  will  appear  at  180“.  The  rear  lobe  usually  is  slightly  narrower 
than  the  front  or  major  lobe  because  the  back  plate  is  laurger  than  the 
array  or  diapiuragm  designed  to  radiate  in  the  forward  direction.  If  one 
assumes  conservatively  that  the  reaur  lobe  has  the  same  beam  width  as  the 
major  lobe,  then  a  correction  to  the  directivity  index  from  Eq.  (2)  is 
easily  estimated  from  the  number  of  decibels  that  the  reau:  lobe  is  below 
the  major  lobe.  Such  corrections  are  shown  in  Table  2.  Minor  lobes  to 
the  rear  at  angles  other  tham  130®  usually  are  small  enough  to  neglect. 


Table  2.  Correction  (in  decibels)  to 
calculated  directivity  index  as  a 
function  of  rear  lobe  height. 


Rear  lobe  down 

(dB) 

D .  correction 
1 

(dB) 

10.0 

-0.4 

12.5 

-0.3 

15.0 

-0.2 

17.5 

-0.1 

20.0 

<-0.1 

Radiation  Resistance 

The  relationship  between  the  radiation  resistance  R  and  the  directiv¬ 
ity  factor  R„  of  a  baffled  piston  can  be  useful  in  ascertaining  the  limit 
of  validity  of  Eq.  (2) . 

The  diffraction  constant  D  of  any  treuisducer  is  given  by  [19] 

=  ERgC^r/Ck^pc)],  (5) 

vdiere  R  is  the  radiation  resistance  in  acoustical  ohms.  Equation  (5)  can 
be  written 

Rg  =  (dVx2)  (pc/R) .  (6) 

For  a  piston  in  a  rigid  baffle,  D  =  2,  and  when  the  piston  is  large, 

R  =  pc/A  and  the  acoustic  load  becomes  largely  resistive.  Then  Eq.  (6) 
becomes , 

R-  =  4ttA/x2,  (7) 

y 

or  the  equivalent  of  Eq.  (2).  Thus,  Eqs.  (2)  and  (7)  are  valid  approxi- 
cations  whenever  a  baffled  transducer  has  a  specific  acoustic  impedaince 
load  of  pc. 
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Equation  (6)  can  also  be  used  for  other  transducer  configtixctions, 
where  R  is  available  from  a  number  of  books  and  D  is  available  from 
Henriquez's  paper  [20]- 


Line  Sources 

The  theoretical  values  of  the  directivity  index  of  line  sources  have 
been  well  known  from  the  work  of  Stenzel  II]  and  Molloy  12] .  However,  the 
analysis  or  calculations  of  botli  authors  have  been  marraU  by  errors. 

Stenzel  [17]  has  pointed  out  some  errors  in  Molloy *s  paper  12].  Stenzel, 
in  turn,  has  erred  in  his  original  work,  and  both  a  second  edition  by 
Brosze  [12]  cmd  translations  into  English  [10,11]  contain  additional  errors. 
Noiie  of  these  latter  errors  are  significant,  but  a  resume  and  correction 
seem  in  order. 

Thirty  yeaurs  ago,  the  terminology  was  not  consistent  between  "direc¬ 
tivity"  euid  "radiation,"  and  between  "factor"  and  "index."  Also,  by  what¬ 
ever  name,  the  directivity  factor  was  the  reciprocal  of  the  modern  paura- 
meter.  Consequently,  in  what  follows,  the  current  definition  of  directivity 
factor  and  appropriate  inversions  will  be  used. 

The  exact  expression  for  the  directivity  factor  Rg  of  a  line  source  is 
given  in  all  sources  as  equivalent  to 


sin2(kL/2) 


2 


kL 

f  sin  t 


(kL/2)2 


dt. 


(9) 


vdiere  k  is  the  wave  number  amd  L  the  length  of  the  line.  Stenzel  [1]  de¬ 
rived  Eq.  (8)  from  the  case  of  a  rectangular  piston  in  an  infinite  rigid 
baffle  (but  radiating  on  both  sides  of  the  baffle)  letting  the  width  uf 
the  rectamgle  approach  zero. 

Stenzel  evaluated  Eq.  (8)  numerically  by  using  tables  for  the  sine 
integral,  amd  analytically  by  using  the  approximation  for  a  sine  integral 


X 


0 


sin  t 

- dt 

t 


n 

2 


cos  X 


sin  X 


X 


(9) 


Substituting  this  r^-proximation  in  Eq.  (8) ,  using  trigometric  identities 
for  doable  angles,  and  reaurranging,  produces 


2  2  sin  kL  4  cos  kL 

kL  kL  (kL)3 


(10) 
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12 


It  should  be  noted  that  the  vertical  scale  in  Fig.  5  is  half  that  of 
Figs,  3  euid  4,  so  as  to  show  the  differences  among  the  various  approxi¬ 
mations  . 

Clearly  the  fovir  terms  are  a  poorer  approximation  than  the  two  terms, 
and  the  two  trigometric  terms  amd  the  errors  in  them  can  be  neglected. 

The  one-tem  approximation, 

Rq  =  KL/it  =  2L/X  (11) 

is  commonly  used  in  transducer  analysis,  but  at  low  frequencies  or  for 
short  lines  the  two  terms. 


• 

ir 

% 

2 

-1 

2L 

X 

1  .  -  -  - 

kL 

{kL)2 

X 

ii2l. 

should  be  used. 

In  terns  of  tiie  directivity  index  =  10  log  R^,  the  errors  cire 

<0.5  dB  for  L  >  0.8X  1  ^ 

V  for  Rg  =  2L/X 

<0.2  dB  for  L  >  2. OX  J 

<0.2  dB  for  L  >  0.5X  for  R»  =  (2L/X)  (1  -  X/Tr“L)"^ 

o 

The  line  transducer  has  no  baffle  conditions  or  rear  lobes  to  be  con¬ 
cerned  with.  The  length  L  is  determined  from  baaro  widths  exactly  as  the 
length  of  a  rectangular  piston;  or  from  the  known  physical  length  that 
includes  a  half- interstitial  space  at  each  end. 


Siiaded  Transducers 

All  of  the  foregoing  approximation  theory  is  based  on  the  cases  of  j 

pistons  or  lines  that  have  uniform  response  over  the  entire  etrea  or  length.  j 

Many  sonar  transducers  are  shaded.  Thrt  is,  the  vibration  an^litudc,  when  I 

transmitting,  is  a  maximinn  at  the  center  and  tapers  off  to  some  lower  val-  | 

ues  toward  the  periphery  or  end.  The  purpose  is  to  suppress  the  side  lobes,  j 

but  an  associated  effect  is  to  widen  the  major  lobe.  Cam  the  approximations  I 

for  uniform  radiators  be  applied  to  shaded  trauisducers?  1 

1 

The  effect  of  suppressing  the  minor  lobes  can  be  estimated  from  Table  1.  J 

In  the  most  extreme  case  of  a  piston  pattern  with  no  minor  lobes,  the  maxi-  J 

mum  correction  is  16%,  or  a  0.6-dB  increase  in  the  directivity  index  of  a 
corresponding  uniform  piston.  For  a  shaded  line,  the  maximum  correction 
is  a  0.3-dB  addition.  Other  corrections  can  be  estimated  within  very  small 
errors . 

Hie  widening  of  the  major  lobe  reduces  to  the  question  of  whether  the 
relative  shape  of  the  lobe  remains  the  same,  or  whether  the  major  lobe 
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approximates  that  radiated  from  a  smaller  uniform  source.  In  addressing 
these  questions,  two  types  of  shading  functions  were  investigated — where 
shading  function  pertains  to  the  mathematical  description  of  the  sensi¬ 
tivity  of  a  radiator  as  a  function  of  the  distance  from  the  geometric 
center. 

The  first  was  a  "linear  taper,"  where  che  sensitivity  varies  linearly 
from  the  majiimum  at  the  center  to  zero  at  the  end  or  periphery.  The  pat¬ 
tern  of  such  a  line  is  given  by 


fsin  [(KL/2X)sin 


Pe  = 


(itL/2X)sin  0 


(13) 


The  first  minor  lobe  of  this  pattern  is  26.6  dB  down. 


The  second  function  is  that  for  a  "binomial  line."  This  function  is 
used  with  a  line  array  of  point  elements  whose  sensitivities  are  propor¬ 
tional  to  the  coefficients  in  the  expansion  of  a  bincmiial  runction 

(x  +  y)^  where  n  is  the  number  of  elements.  If  the  element  spacing  is 
a  half-wavelength,  the  shading  is  perfect  or  there  are  no  minor  lobes. 

For  the  investigation  here,  n  was  10.  The  optimum  half-wavelength  spacing 
was  used,  thereby  simulating  a  five-wavelength  continuous  line.  The  co¬ 
efficients  in  the  expansion  of  (x  +  y)®  are  1,  9,  36,  84,  126,  126,  84, 

36,  9,  and  1.  When  normalized  so  that  their  sum  is  one,  these  coefficients 
become  0.002,  0.018,  0.070,  0.164,  0.246,  0.246,  0.164,  0.070,  0.018,  and 
0.002.  The  pattern  of  such  a  line  array  with  half-wavelength  spacing  is 
given  by 


p  =  2[o.246  cos(0.5tt 
+  0.070  cos (2.511 
+  0.002  cos  (4.511 
The  I'elative  beam  widths  of  the 

3  dB 
6  dB 
10  dB 


sin  9)  +  0.164  co5(1.5ir  sin  9) 
sin  9)  +  0.018  cos(3.5ii  sin  6) 
s.  n  0)  ]. 

patcem  of  Eq.  (14)  cure 
down  0.70 

down  1.00 

down  1.23 


(14) 


Comparing  these  relative  beam  widths  with  those  of  uniform  radiators,  it 
is  seen  that  the  two  cases  are  similar,  but  not  identical. 

To  obtain  a  quantitative  effect  on  the  directivity  index,  a  uni.'’orm 
line  length  auid  a  linecirly  tapered  line  length  were  chosen  so  that  their 
patterns  had  the  same  6-dB-down  beam  width  as  the  bincxnial  line.  These 
lengths  turned  out  to  be  2.4<^X  and  3.59X,  respectively.  The  major  lobes 
t-f  the  three  patterns  are  s’.ow,.  in  Fig.  6. 

Graphical  calculation  ,f  the  directivity  index  using  the  patterns 
shown  as  the  major  lobe  of  a  llr  ?  so-orce  (having  a  toroidal  pattern) ,  and 
neglecting  all  minor  lobes,  gives  the  following: 
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Uniform  line 


6.6  dB 


Linearly  tapered  line  6.5  dB 


Binomial  line  array 


6.5  dE 


The  theoretical  value  for  the  uniform  line  is  6.9  dB  for  the  whole  pattern 
and  7.2  dB  for  only  the  major  lobe.  The  0.6-dB  error  in  tho  graphical 
calculation  is  not  unexpected  in  light  of  the  data  in  Fig.  2  and  the 
imprecision  in  measuring  levels  at  angles  near  the  pattern  axis. 

If  the  patterns  in  Fig.  6  are  taken  as  those  for  a  circular  piston, 
the  graphical  calculations  yields  the  following  directivity  i.‘dexes: 

Uniform  line  19.4  dB 

Linearly  tapered  line  19.6  dB 

Binomial  line  aurray  19.5  dB 

The  conclusion  here  is  that  the  differences  in  ■^Le  patterns  of  uniform 

and  shaded  radiators  are  less  them  measurement  error  insofar  as  the  effect 
on  the  directivity  index  is  concerned.  Consequently,  the  technique  of 
measuring  the  6-dB>down  beaun  width  and  then  calculating  the  directivity 
index  as  if  the  transducer  were  uniform  is  acceptably  accurate  for  shaded 
transducers. 


unifors 


Fig.  6.  Major  lobes  of  patterns  for  a 
uniform  line,  a  lineairly  tapered  line, 
amd  10-eleroent  binomially  shaded  line 
array. 


Arrays  with  Mutual  Coupling 

One  other  case  of  nonuniform  radiators  should  be  mentioned.  In  large 
low-frequency  aurrays  where  the  elssent  spacing  is  a  small  fraction  of  a 


wavelength,  there  is  mutual  interaction  or  coupling  .unong  the  elements. 

That  is,  the  radiation  impedarce  of  one  element  is  affected  by  the  vibra¬ 
tion  of  neighboring  elements.  The  elements  usually  are  vibrating  at  their 
resonance  frequency  and  thus  are  sensitive  to  ciny  effect  on  their  radiation 
impedance.  The  result  can  be  nonuniform  vibration  siuong  the  elements  due 
to  this  mutual  coupling,  which  would  invaJidate  the  directivity  index 
computation  based  on  uniform  pistons  or  ?ines.  As  in  the  case  of  shaded 
transducers,  it  becomes  a  question  of  whither  the  measured  pattern  is  simi- 
Icu:  enc-ugh  to  that  of  an  equivalent  uniiurm  eurray.  Unlike  shading,  mutual 
coupling  it  an  tmintentional  2uid  undesirable  effect.  When  a  pattern  is 
significcntly  affected  by  mutual  couplir»g,  the  problem  usually  is  that  of 
correcting  the  cause  rathe-  than  measuring  L.he  result. 


Conclusion 

From  both  experiment  and  theory,  it  is  apparent  that  the  directivity 
index  of  any  ordinary  transducer  can  be  obtained  from  calculations  based 
on  known  configuration  and  dimensions  or  beam-width  measurements  with  a 
degree  of  reliability  and  acciuracy  that  is  no  worse  than  any  measurement 
technique.  A  censervat-ive  limit  for  the  validity  of  such  calculations  is 
that  the  minimum  trans-*  -r  dimension  be  one  wavelength. 


Ackn(h.<cugement 

The  author  is  indebted  to  James  D.  George  for  computer  calculations 
used  in  this  study. 


References 

{1]  H.  Stenzel,  Leitfaden  zur  Berechnung  von  Schallvorgangen  (Julius 
Springer,  Berlin,  1939);  also  republished  by  J.  W.  Edwards  in  1944 
in  the  United  States  under  the  authority  of  Alien  Property  Custo¬ 
dian  License  Ito.  A-491. 

[2]  C.  T.  Molloy,  "Calcnilation  of  the  Directivity  Index  for  Vauriovrs  Types 
of  Radiators,"  J.  Acoust.  Soc.  Am.  20,  387-405  (1948). 

[3]  P.  H.  Kendig  amd  R.  E.  Mueser,  "A  Simplified  Method  for  Determining 
Transducer  Directivity  Index,"  J.  Acoust.  Soc.  Am.  19,  691-694  (1947). 

[4]  C.  E.  Green  amd  J.  R.  Roshon,  "Directivity  Factor  Computer  for  Elec¬ 
troacoustic  Tramsducers,"  MEL  Report  1196,  U.  S.  Navy  Electronics 
Ladx>ratory,  Sam  Diego,  California,  13  Sep  1963. 


16 


Scientific-Atlanta,  Inc.  is  the  manufacturer  of  equipment  for  inte¬ 
grating  the  patterns  of  both  sonar  transducers  and  antennas. 

A.  M.  Young,  "Digital  System  for  the  Measurement  of  Directivity 
Index,"  NRL  Report  7585,  20  Apr  1973  [AD-758  639]. 

R.  F.  Green,  ".Measuring  the  Directivity  Index  of  Underwater  Sound 
Projectors,"  IEEE  Trans,  on  Audio  and  Electroacoustics  AU-21 ,  4C7-412 
(1973) . 

R.  J.  Bobber,  Undemater  Electroacoustic  Measurements  (Naval  Research 
Laboratory,  U.  S.  Government  Printing  Office,  Washington,  D. 

1970),  pp.  83-90. 

Special  slide  rules  for  calculating  sonar  transducer  parameters  have 
been  pror’.uced  by  the  Edo  Corporation,  Sperry-Rand ,  and  The  Raytheon 
Corp. 

H.  Stenzel,  Guide  for  the  Calculation  of  Sound  Processes,  navships 
250-940;  translation  of  reference  1  by  C.  E.  Mongan,  Jr.,  Hay  1947. 

H.  Stenzel,  Handbook  for  the  Calculation  of  Sound  Propagation 
Phenomena,  NRL  Translation  130  by  A.  R.  Stickley,  Nov  1947  [trans¬ 
lation  of  reference  1  into  English). 

H.  Stenzel  and  O.  Brosze,  Leitfaden  zur  Berechnung  von  Schallvorgangen, 
Zweite  Auflage  (Springer-Verlag ,  Berlin,  1958)  [a  second  edition  of 
reference  1,  in  German,  by  Brosze  after  Stenzel ’s  death). 

Reference  8,  pp.  89-90. 

B.  G.  Watters,  "A  Reverberant  Tank  for  Underwater  Measurements," 

J.  Acoust.  Soc.  Am.  53,  357(A)  (1973)  [Paper  II4,  84th  iMeeting 
of  the  Acoustical  Society  of  America,  Nov  1972). 

A.  E.  Reznikov  amd  A.  Ya.  Snytko,  "Problem  of  Measuring  the  Axial 
Concentration  Coefficient  in  Ultrasonic  Radiators,"  Measurement 
Tecnniques  [the  Soviet  Journal  Izmeritel 'naga  Tekhnika  in  English 
translation).  No.  7,  July  1965,  pp.  654-657. 

L.  L.  Beranek,  Acoustics  (McGraw-Hill  Book  Co.,  New  York,  1954), 
p.  112. 

H.  Stenzel,  "Remarks  on  a  Paper  entitled  'Calculation  of  the  Direc¬ 
tivity  Index  for  Vaurious  Types  of  Radiators,’"  J.  Acoust.  Soc.  Am.  24, 
417-41fJ  (1952). 

R.  J.  Bobber,  "The  Effects  of  Element  Packing  on  the  Con5>lete  Radia¬ 
tion  Pattexms  of  Arrays,"  NRL  Memorandum  Report  2206,  18  Jan  1971 
[AD-718  312). 

R.  J.  Bobber,  "Diffraction  Constants  of  Transducers,"  J.  Acoust. 

Soc.  Am.  37,  591-595  (1965). 

T.  A.  Henriquez,  "Diffraction  Constamts  of  Acoustic  Tramsducers , " 

J.  Acoust.  Soc.  Am.  36,  267-269  (1964). 

American  National  Standaurd,  "Procedures  for  Calibration  of  Under¬ 
water  Electroacoustic  Transducers,"  ANSI  SI. 20-1972,  p.  27. 


17 


